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Abstract: In this work we discuss our consistent implementation of the minimal 
model of Universal Extra Dimensions in CalcHEP. We pay special attention to the 
gauge invariance issues that arise due to the incorporation of 5D quantum corrections. 
After validating the implementation we perform a complete study of the tri-lepton 
signature, including a realistic estimate of the backgrounds, for the present LHC 
energy and luminosity. We also derive the expected LHC discovery reach for different 
luminosities, both at ^/s = 7 TeV and 8 TeV. 
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1 Introduction 

Many beyond the Standard Model (SM) theories rely on the existence of extra spatial 
dimensions. For instance, the best candidate to unify gravity and gauge interactions, 
string theory [1], is only consistent in a ten dimensional space-time. On the other 
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hand, the presence of extra dimensions could also have an impact on scales much 
lower than the Planck mass [2] (and this can be consistent with string theory [3]). In 
particular, extra dimensions provide a new perspective from which to address issues 
such as the nature of electroweak symmetry breaking, the origin of the number of 
fermion generations [4] and their mass hierarchy [5], and the supersymmetry breaking 
mechanism [6]. In any case, the underlying physics will depend on several facts 
such as the size and the shape of the extra dimensions, the existence of space-time 
subspaces in which the different fields propagate, the presence of background fields, 
etc. Let us concentrate on the size of these extra dimensions. Although present data 
constraint them, they could be as large as C(TeV)~^. In this case they would be 
accessible at the Large Hadron Collider (LHC), which is already probing this scale. 
Furthermore, since this is also the typical mass of a weak interacting massive particle 
(WIMP), such models could be relevant for dark matter purposes. 

The simplest of these models is the Universal Extra Dimension (UED) scenario, 
which was proposed a decade ago by Appelquist, Cheng and Dobrescu [7]. In this 
model, the geometry consists of the ordinary flat 4D space-time and an additional 
fifth dimension compactified on a circle of radius R. In its minimal version (MUED), 
the SM particles are promoted to 5D fields which propagate in the whole space- 
time. Notice that since gauge couplings are not dimensionless, the model is non- 
renormalizable. Therefore it should be regarded as an effective theory valid up to 
some cut-off A, above which new degrees of freedom that complete the theory have 
to be incorporated. Since the extra dimension is compact, for small values of the 
radius R the theory may be reduced to an effective four- dimensional theory whose 
light spectrum contains the SM. The precise SM spectrum is obtained by imposing a 
Z2 symmetry that breaks the translational invariance along the fifth dimension and 
projects out half of the 4D fermions, allowing for chirality. This symmetry does not 
break Kaluza-Klein (KK) parity, a discrete symmetry that ensures the stability of 
the Lightest KK Particle (LKP), the potential dark matter candidate of the model. 

The spectrum is organized in towers of KK modes that get a common mass 
contribution proportional to their KK number. The degeneracy of these KK masses 
is lifted by the breaking of the electroweak symmetry and by 5D quantum correc- 
tions. In fact, a detailed evaluation of the 4D spectrum is crucial to explore the 
phenomenology of the model. From one side, it determines the sequence and the 
kinematics of the KK chain decays that will be produced at the LHC. On the other 
side, it enters in the evaluation of the LKP cosmological abundance [8], in which 
coannihilation processes are expected to play an important role. 

It is then important to deal simultaneously with the breaking of both gauge and 
space-time symmetries and to understand their interplay with quantum corrections. 
This question has been addressed in several works. In particular, analytic expres- 
sions for the mass spectrum without including radiative corrections were obtained 
in ref. [9]. The effect of 5D quantum corrections was first evaluated in ref. [10] in 
an approximation such that the effects which involve both electro-weak symmetry 
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breaking and radiative corrections were ignored. Since these terms are not gauge 
invariant, these approach is not fully satisfactory. A further step along this lines 
was made in ref. [11], where the 5D corrections (both bulk and brane) were incorpo- 
rated as a series of wave function renormalization factors in the 4D expansion before 
considering the electroweak symmetry breaking. The consistency of this procedure 
has to be analyzed since some of these terms cannot be derived from a local 5D 
Lagrangian. 

The aim of this work is to provide a consistent description of MUED, implement 
it in CalcHEP (a software package [12] for calculating cross-sections and generating 
events) and explore the signatures of this model to study its discovery potential at 
the LHC. 

We have documented our implementation, collecting all the expressions (com- 
plete Lagrangian terms, masses and mixing values after EWSB, etc.) and provided 
detailed tables describing the notation used for all particles and parameters. We have 
also made the model publicly available at the High Energy Model DataBase (HEP- 
MDB) in both Feynman-'t Hooft and unitary gauges at https:/ /hepmdb. soton.ac.uk/ 
under the name "MUED-BBMP" and have used the power of HEPMDB to perform 
scans over the MUED parameter space. Using this setup we have studied the LHC 
potential to explore the MUED parameter space using a tri-lepton signature, which 
we found to be one of the most promising signatures sensitive to MUED. There are 
several reasons for this. One is related to the compressed spectrum of MUED which 
leads to soft particles in the final state, among which the leptons are the most con- 
venient objects to be identified and dealt with. Another reason for the advantage of 
the leptonic signatures, including the tri-lepton one, is the specific mass hierarchy 
of the KK particles in MUED which provides a large branching ratio of leptonic 
decays from initial KK states; we demonstrate this below. Finally, the tri-lepton 
signature has quite a high signal-to-background ratio. In addition, we have identi- 
fied optimal kinematical cuts that further enhance the statistical significance of the 
signal. The results obtained are discussed in this article and compared with the pre- 
vious studies [13-17] devoted to the MUED collider phenomenology. The analysis is 
done in the region compatible with present indirect constraints (electroweak preci- 
sion measurements [18], Higgs boson production and decay data [19, 20] and Dark 
Matter [11]). 

The detailed content of the paper is the following: in section 2 we present our 
theoretical framework and provide the detailed MUED Lagrangian including the 5D 
quantum corrections. We also describe the spectrum for two benchmark points. In 
section 3 we deal with the CalcHEP model implementation and its validation. The 
phenomenology of the model is presented in section 4. We study in detail there 
the tri-lepton signature, paying special attention to the background evaluation, and 
make a realistic estimate of the LHC discovery reach. In section 5 we summarize our 
results. In appendix A we give the explicit relations for masses and mixings in the 
neutral bosonic sector, while in appendix B we describe the implementation of four 
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gluon vertices in CalcHEP. 



2 Theoretical framework 



In this section, we start by reviewing the most important features of the MUED 
model [7]. We then explain our formalism in detail and show that it reconciles gauge 
invariance, electroweak symmetry breaking and 5D loop corrections. 

2.1 The MUED model at tree level 

In the MUED Model, all SM particles are promoted to 5D fields which propagate in 
the bulk of an /'L2 orbifold^ with length ttR. Because of the periodicity of S"^, the 
fields can be expanded in Fourier modes exp (^) with n an integer. However, half 
of these states are projected out under the action of the orbifold. This procedure 
allows for fermion chirality. At the end of the day, the y-independent n = modes 
reproduce the 4D SM particle content while in turn each SM particle is associated 
to a tower of Kaluza-Klein (KK) partners. 

2.1.1 Decomposition of the fields 

Under a Z2-transformation, a scalar field (f) can be even or odd, i.e </> {—y) = ±</> [y), 
because its kinetic term contains two derivatives. Therefore its expansion in KK 
modes is either 
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We see that an odd scalar does not have a zero mode and hence is not suitable to 
represent the Higgs field. Furthermore, the zero mode 0'-*^^ of an even scalar field 
comes accompanied by the infinite tower of (^'•"^ states. We note that the factor 
1 / \/ttR accounts for the dimensionality of the field in 5D and ensures that the KK 
fields are canonically normalised in 4D. 

The decomposition of a fermion can be deduced from its kinetic term ip'y^'^dM^', 
which contains ip'-f^d^ip. This term is invariant under Z2 if ip{—y) = ±^^ip{y). If 
the sign is + (resp. — ), then the zero mode of the fermion is right-handed (resp. 
left-handed) . 
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^We denote the extra dimension by y or and the 4D coordinates by . The metric is 
diag(+, —,—,—,—). Under a Z2-transformation, y — > —y. 
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Note 



where the subscript L, R refers to the zero mode chirahty, and Pl,r 
that all the ip^^l^ are Dirac fermions while the zero mode is chiral. Thus, there are 



four chiralities per KK level : V'L^^ ^l^^ V'/Ti ^Rr' Only when necessary will we 
specify these chiralities. Otherwise, it should be understood that ip^j^^j^ with n > 1 
stands for a Dirac fermion. 

Finally, a gauge boson Am, M = (/i,5), couples to a fermion via the term 
iJj'J^^DmiP ^ ipl^'^ Am^P- Consequently, is even while A^ is odd 



^5 {x,y) 




V2 J] COS 



X) sm 



n>l 



n>l 

ny_ 
R 



ny 
R 



(2.4) 



2.1.2 The MUED Lagrangian 

It is simply obtained by embedding the SM Lagrangian into the 5D geometry: 



+i QT^^DmQ + I LT^'DmL + i uT^'Dmu + i dV^'Dud + i eV^'Due 
+H^H^H - Xh [H^H) ^ - (^yuQuH + ydQdH + yMH + h.c.) , (2.5) 

where G, W and B are the gauge bosons of respectively SU(3), SU(2) and U(l); H 
is the Higgs field; Q and L are the quark and lepton doublets; and u, d and e are 
the up-type quarks, down-type quarks and electron-like singlets. All fields should 
have an understood (x'*, y) dependence and there is an implicit sum over fermion 
generations. 

The Dirac matrices are F^^ = (7^, Z7^) and the covariant derivative is 



D 



M 



with Y the hypercharge, a* the Pauli matrices, the Gell-Mann matrices and where 
the 5D gauge couplings are related to the 4D ones by g\ = giVirR. Finally, the 
factors yu, yd and ye are the Yukawa matrices and H = ia^H*. 

It is relatively straightforward to expand the various fields into KK modes and 
integrate the Lagrangian (2.5) over the compact extra dimension in order to study 
the effective 4D theory: = j^^ ^^{y)dy. Due to the kinetic terms along the fifth 
dimension, the KK particles receive a geometrical contribution to their mass 
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"^fermion(n) = ^0 + — (2.6) 

K 

where mo stands for the mass of the corresponding zero mode, that is, the SM mass. 
These expressions are modified when quantum corrections are included. 

2.2 MUED at one loop: mass corrections and gauge invariance 

According to (2.6), the lightest KK modes are the n = 1 photon and gluon, which 
have the same mass. Other degenerate states are the two KK Dirac fermions associ- 
ated with each quark or lepton. In fact, when the compactification scale is much 
larger than the Higgs vev v, the entire spectrum at a given KK level is approximtely 
degenerate. Under these circumstances, one expects radiative corrections to play an 
important role when studying the phenomenology of models with a small compact 
extra dimension^. 

Quantum corrections to the KK masses arise from the breaking of the 5D Lorentz 
invariance. They have been studied [10, 22-24] in detail in the v ^ limit. There 
are two types of quantum corrections to the KK masses: the so-called bulk and brane 
contributions. The first ones are generic in a compactified theory since 5D Lorentz 
invariance is broken at long distances by the compactification. These corrections 
are intrinsically non-local (in 5D), modify only the gauge boson masses and are 
independent of the KK level. 

The second ones [22] have their origin in the breaking of the translational sym- 
metry at the orbifold fixed points. Brane corrections are divergent in the cutoff scale 
A and require brane-localised counter-terms to renormalise the action. In the effec- 
tive 4D theory, these corrections are proportional to the usual Kaluza-Klein mass ^ 
and depend logarithmically on A. 

Summing both bulk and brane contributions, one has 

"'^vector (n) — (^v ' ^ ^ 

A - !L 

'J'^fermion (n) '2/ j-, ; 

Jri 

2 

'^"^scalar (n) ^ ^ ' C^'^) 

for n > 1, where the factors a^, af, as, are dictated by the quantum numbers of 
the particle (see [10] for the complete set of corrections). The first terms represent 
the brane corrections, while the constants stand for the bulk corrections. For 
example, for the quarks doublets, one has 

2 27 2 1 2^ . A2 



= 16^ l'^^ + 16^^ + 16^^ J 

where the gi are the gauge couplings for the different SM gauge groups, A is the UV 
cutoff, and /i is the running scale. Comparatively, the bulk corrections for the U(l) 



^The importance of quantum corrections in this context was first pointed out by Hosotani [2f]. 
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jauge boson B are 

^ 39 C (3)^1^ 1 
2 167r4 i?2 



As mentionned in the Introduction, adding these corrected mass terms to the 
Lagrangian (2.5) exphcitly breaks gauge invariance in 5D. In the remainder of this 
section, we model the brane corrections with wave-function normalisation factors as 
described in [11], and proceed with the compactification. The strategy is to focus on 
the quadratic part of the Lagrangian in order to identify the mass eigenstates, the 
Goldstone bosons and Fadeev- Popov ghosts. Note that the quadratic part of the 5D 
Lagrangian is diagonal in the Kaluza-Klein basis, see e.g [9]. We comment on the 
bulk corrections at the end of this section. 

2.2.1 The case of a scalar field 

Before turning to the full MUED Lagrangian at one loop, let us consider a scalar 
field $ in 5D. The key point of our setup is to recognize that the brane corrections 

2 

are proportional to the KK mass ^ and hence can be incorporated into the kinetic 
term along the fifth dimension. We introduce normalisation functions by splitting 
the 5D part from the 4D part in the Lagrangian 

-^^ = I -l^^ l^5$|' + ml |$|' . (2.8) 

Following section 2.1, we decompose the field in KK modes and integrate over the 
fifth dimension, ^4 = j^^ (y) dy. In 4D, the spectrum consists in a tower of KK 
scalar fields 0'-"^ with masses = ml + Z,^^. It is straightforward to identify 
Zq, = 1 + as so that the total mass matches the tree-level mass (2.6) plus the brane 
correction (2.7). 

Let us comment on this procedure. First of all, the Lagrangian (2.8) is not 
Lorentz covariant. However, since the orbifold breaks translational invariance in 5D, 
we emphasize that this Lagrangian is phenomenologically relevant as far as brane 
corrections go^. Also, in principle, loop diagrams of the "brane" type yield KK- 
number violating mass mixings since is not conserved in these processes (see [10], 
eq. (34)). This means that in the 4D effective theory, the KK states are no longer the 
mass eigenstates. One should then diagonalise the infinite mass matrix. However, 
the off-diagonal KK-mixing terms in this matrix are small compared to the diagonal 
mass terms and so, to a very good approximation, we can disregard them and retain 
only KK-number preserving mass corrections. In 5D, this approximation manifests 
itself through the fact that our Z^ function is independent on the extra dimension, 
and hence KK number is conserved. In other words, instead of considering brane- 
localised counter-terms, we have effectively spread them throughout the compact 
dimension. 

^Alternatively, this is equivalent to rescaling R R/ \/ Z^. 
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2.2.2 The MUED Lagrangian at one loop 

The same idea of sphtting the fifth direction from the Minkowski 4D can be appUed 
to the whole Lagrangian (2.5). Whenever there is a contraction of Lorentz indices, 
that is for all the kinetic terms, we introduce a Z function. 

This procedure preserves gauge invariance. For instance, consider a field strength 
jp<^MN -pa^^^ It is gauge invariant in 5D by construction. Due to its antisymmetry, it 
reduces to F'^^^F"'^^ + 2F^^F"-^^ . Each of these two terms is clearly gauge invariant 
on its own, therefore splitting them explicitly in the Lagrangian and adding a factor 
Z to the second one does not affect gauge invariance in 5D. The same applies to the 
Higgs and the fermions kinetic terms. 

We are now ready to write our 5D MUED Lagrangian and analyse it: 

Zc^a r^afj. , Zb DM i 

^*-^m5^ 5 + ^-^m5-D 5 + 

-ZqQ-i'^D^Q - ZlL^'D.L - Z^u^'D.u - Z^d^^D.d - Z.e^^'D.e 
+H'^H^H - Xh [H^nf - (y^uH + y^dH + y^H + h.c.) . (2.9) 

We will focus on the quadratic part of the Lagrangian in order to identify the mass 
eigenstates of the 4D theory. Also note that the quadratic part of the Lagrangian is 
diagonal in the KK indices. Therefore, the strategy is the same as for MUED at tree 
level, see [9]: we KK expand the fields and perform the compactification. We then 
diagonalise the quadratic part of the Lagrangian in 4D for a given level n. After 
identifying the Goldstone bosons, we perform the gauge fixing and incorporate the 
Fadeev-Popov Lagrangian. 

The bosonic sector 

For the reader interested in the details of the computation, an example of diag- 
onalisation of the Lagrangian (2.9) is presented in appendix A where we consider 
the subgroup U(l)xU(l) spontaneously broken by the Higgs. As shown therein, 
the identification of the Goldstone bosons is tedious, though not conceptually dif- 
ficult, and the gauge fixing proceeds as usual. We list here the key results of the 
diagonalisation procedure. 

Once the Higgs field acquires a vev, the bosons B and mix. However, the 
KK states and ^^^(n) 

mix with an angle that depends on n, Zb and Z^. Con- 
sequently, the n > 1 mass eigenstates cannot be identified with a KK photon or a 
KK Z-boson. We call the mass eigenstates PjT^ and V^T^ and their masses are 

1 Tl^ 

^{p,v}in) = 9 T ^91)+^ (Zb + Zw) 




91) + ^,{Zw-Zb)) +-glglv' 



. (2.10) 
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Likewise, the fifth components B^"'\ W^^""^ and the Higgs pseudo-scalar also 

(n) 

mix. Two of the mass eigenstates are the Goldstone bosons associated with and 
Y^"'\ The third scalar mass eigenstate is a Higgs neutral field Oq"'' which mass is 



i?2 ' A\Zb Zw 



(2.11) 



On the other hand, the gauge bosons Wlj and Wl,^ give rise to a tower of KK 
states PV/t^"^ since they are orthogonal to the breaking of SU(2). Their mass is 
simply 



n 



2 



m^(„) = + Zw^, ■ (2.12) 

The associated Goldstone boson G^''"^ is a combination of their fifth component 
W^^'^'^ with the Higgs pseudo-scalars x^''"^ ^"^^ The other linear combination 

is a physical charged Higgs scalar a|^^ whose mass is 

ml{n) = §^ "^VF(n) • (2-13) 

The only remaining degree of freedom of the Higgs, which we call H^'^\ acquires a 
mass 

^l(n) = M| + Z^— . (2.14) 

Its zero mode corresponds to the SM Higgs boson. 

Finally, the gluon has a KK tower G^^"'' with mass 



2 



I^Gin) = ^G-^ , (2.15) 

and its associated Goldstone bosons are the fifth components 

Therefore, at the n*^ KK level, there are five gauge bosons: G^^\ P/i"'', v/r\ 
W/t^"'^', and four scalars: the Higgs H'^'^^ and the three scalars Oq"^ and a^]!\ We 
point out that among the scalars, only the Higgs field possesses a zero mode. 

Once we have identified the Goldstone bosons of the theory, the gauge-fixing 
Lagrangian is determined very easily, see for example the end of appendix A. The 
next step is to include the ghosts fields. As an illustrative example, let us focus on 
the SU(3) sector. The 4D gauge-fixing Lagrangian 



2e 

can be promoted to a 5D Lagrangian 

^GF = {d^'G^, - iZcd^GlY = -y^F^F^ , (2.16) 

where we have introduced the gauge-fixing functions Fq. Because this Lagrangian 
is quadratic, it cannot give rise to any KK-number mixing. The next step is to 
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introduce the 5D ghosts fields and c", which are Z,2-even. The Fadeev-Popov 
Lagrangian is 

where a° is the parameter of the gauge transformation. We find 

For the SU(2)xU(l) sector of the theory, particular care has to be taken since the 
symmetry is spontaneously broken. We emphasize however that the Fadeev-Popov 
procedure in this case does not bring any additional complication compared to the 
SM. 

The fermionic sector 

We now proceed to the diagonalisation of the fermionic sector of the theory. The 
left and right fermions do not receive the same corrections since they are respectively 
doublets and singlets of SU(2). Recall that each fermion of the SM has two associated 



infinite towers of KK Dirac fermions if)^^'^ and i^^lf' ■ 
breaking, they mix with the following matrix 



' '7 n 
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After electroweak symmetry 



(2.17) 



where is the zero-mode mass and we labelled the left and right radiative correc- 
tions as Zl At the n**^ KK level, the mixing is 



with 



tan (2a^(„)) 



-7-" cos q;^(„) smQ!^(„) 
7^sina^(.„) cosa^(„), 

{ZL + Zn)n/R ' 



for n > 1 



and where i/jI"^^ are the mass eigenstates after EWSB, with masses 
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(2.18) 



(2.19) 



(2.20) 



The same applies to all fermions except for the neutrinos. 

Let us briefiy comment on the corrections to the top quark mass. Because 
the Yukawa coupling is large, the top receives an additional correction through the 
Higgs, see [10]. These corrections can be included as extra contributions to the 
wave-function normalisation. However, because of gauge invariance, the left-handed 
bottom quark also receives a Higgs-induced correction. The corresponding terms in 
the Lagrangian are 



Ztl T^'D.T - ZtRh'D^t - Zd b^'D.b 
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where is the common correction for right-handed quarks, and Z^l and Zt^ contain 
the new contributions. Consequently, the bottom masses depend on Ztl instead of 
the naive expectation [Zq). 

In summary, at the n*^ KK level, there are two Dirac fermions labelled by an 
index (1 or 2). Their masses are given by (2.20) and depend on the left and right 
corrections. 

2.2.3 Comments on the bulk corrections 

As mentionned at the beginning of this section, the bulk mass corrections arise when 
one of the internal lines of the 5D loop diagram winds around the extra dimension. 
They affect only the gauge boson masses and are independent of the KK number, that 
is, they shift the entire tower of KK states except for the zero mode. Interestingly, the 
bulk corrections computed in [10, 23] seem to be independent of the regular isat ion 
scheme [25-28]. For instance, in latticised (or deconstructed) versions of 5D theories 
[29-31], which are often considered UV complete, the same results were found. 

Because they do not depend on the KK number, the bulk corrections cannot 
be incorporated within a wave-normalisation function in the manner of the previous 
paragraphs. As it turns out^, they can be modeled by non-local operators (Wilson 
lines) in 5D. The reason is that since one internal line winds around the fifth di- 
mension, such loops cannot be shrunk to a point and hence they are non-local from 
the 5D point of view (but local in the effective 4D theory). Even though we believe 
that this issue deserves further attention, it goes beyond the scope of this work to 
consider non-local operators. 

Quantitatively, our assumption is justified because the bulk corrections are neg- 
ligible. The results derived in [10] are 



5m 



2 ^ 39 C(3)(7f 1 
2 167r4 B? 



dm^in) - 2 16vr4 forn ^ . 



3C(3)^7| 1 



2 167r4 i?2 

Numerically, for all n > 1 and for phenomenologically acceptable values of i? (~ 1 
TeV), the bulk corrections are very small compared to the brane ones with the only 
exception of the B'^^'' boson. However, because gi is so small, it turns out that both 
brane and bulk corrections are in any case negligible compared to the tree-level mass 
l/i?2 for B^^\ Therefore, in our implementation and phenomenological studies, we 
have disregarded the bulk corrections. 

2.3 Detailed analysis of the spectrum 

Armed with the masses computed in section 2.2, we are ready to plot the spectrum 
and analyse its features. For each 5D field, the Z function is matched with the 



*We warmly thank A. Pomarol for very valuable discussion on this matter. 



- 11 - 



results obtained in [10]. In order to efficiently compare our results with the existing 
literature, we fix the cutoff A = 20R~^ and the running scale n = . Furthermore, 
we choose two benchmark values for the compactification scale = 800 GeV and 
R^^ = 1500 GeV. As we will explain in section 4, these points are close to the edge 
of the allowed region. The spectrum of the first KK level is shown in Figs. 1 and 2 
for a Higgs mass niH = 120 GeV and these two values of R~^. 
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Figure 1: The first KK level of the MUED spectrum for R-"" = 800 GeV, ulh = 120 GeV, 
A.R = 20 and = 1, at tree level (left) and one loop (right). 
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Figure 2: The first KK level of the MUED spectrum for i?"^ = 1500 GeV, mu = 120 
GeV, KR = 20 and ^iR = 1, at tree level (left) and one loop (right). 



First of all, the degeneracy of the tree-level spectrum is significantly lifted once 
we include corrections. As expected, the lightest KK particle (LKP) is the gauge 
boson P*^^) which is at 99% composed of the B^^^ boson because the electroweak 
scale is much smaller than the radius. It is stable and hence is a good dark matter 
candidate. The next-to-lightest KK particle (NLKP) for the parameters of Fig. 1 is 
a charged KK lepton, while as seen in Fig. 2, the KK lepton and the scalar are 
almost degenerate. Also note that the V^*^^-* boson consists essentially of the VT^^^) 
boson and hence the three SU(2) bosons are still degenerate at one loop. 

The heaviest KK state is the gluon because 5^3 is large and the brane corrections 
are positive. This is also the reason why the quark masses are lifted by a greater 
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amount than the leptons. The LHC being a hadron colhder, as KK quarks and 
gluons are produced we expect them to decay into the several lighter KK species 
plus SM particles. Notice however that the corrections in the top quark masses by 
the Higgs are negative and proportional to the top Yukawa coupling thus pushing 
these masses down. Remarkably, for = 1500 GeV there is just enough gap for 
the second top tg^-* to decay into a SM top plus P^^^ while at R^^ = 800 GeV this 
decay is forbidden. 



1.3 

^1.2 

1.1 

1.0 
( 

AR 

Figure 3: Evolution of the n = 1 masses with the cutoff AR at fixed R'^ = 800 GeV and 
fiR = 1. 

It is interesting to study the dependence of the n = 1 spectrum on the external 
parameters of the model, namely the cutoff A and the radius R. In Fig. 3, we show 
the logarithmic evolution of the rescaled masses M(A) x R with the cutoff AR for 
a fixed R^^ = 800 GeV and fiR = 1. Because the corrections scale with the gauge 
couplings the gluon mass varies over 20% from AR = 5 to 40 while the P mass is 
almost independent on A. Notice that there are some crossings in the spectrum. 
This means that depending on which value we choose for the cutoff, some decay 
channels may appear that do not exist otherwise. Let us emphasize again that the 
cutoff dependence cannot be subtracted from the theory since the MUED model is 
non-renormalizable. Therefore, there is no theoretical motivation for choosing one 
value of A over another. 

In Fig. 4, we plot the rescaled masses M{R) xR a.s functions of R for a fixed ratio 
A//i = 20. Because we incorporated the quantum corrections in a gauge invariant 
way, we see that the rescaled masses exhibit some dependence on R, whereas there 
wouldn't be any had we added the corrections by hand. Most notably it is 
the electroweak sector of the theory, the Higgs, the W bosons, the tops and the 
scalars and a± which acquire the strongest dependence. When R~^ is close to the 
electroweak scale, the interplay between loop corrections and symmetry breaking is 
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Figure 4: Evolution of the n = 1 masses with the compactification scale at fixed 
K/^l = 20. 



considerable and leads to some crossing in the spectrum. For much larger values of 
i?"^, both scales decouple, the spectrum flattens and we recover the behaviour found 
in [10]. 



3 Model implementation 
3.1 Model file 

We have implemented the MUED model presented in section 2.2 into the LanHEP 
package for automatic evaluation of Feynman rules [32]. Its output files are com- 
patible with the CalcHEP/CompHEP matrix element generator [12]. To date, there 
exist two public implementations of the MUED model corrected at one loop but 
none of them contains a KK Higgs sector. The first one is realised as a FeynRules 
model file [33], while the other consists in a CalcHEP input file created manually 
in the unitary gauge [34]. The LanHEP package was updated in [35] in order to 
treat extra-dimensional models. Consequently, we were able to implement the model 
directly in 5D and have LanHEP integrate out the extra dimension automatically. 
Our model file contains three main parts: the declaration of parameters, variables 
and fields, the reconstruction of the 5D fields in terms of the mass eigenstates and 
finally the 5D Lagrangian. 

The first step is to declare the external variables of our model, which we list in 
Table 3.1. Except for the compactification scale -R~^, the cutoff A and the running 
scale /i, the rest of the parameters are those of the SM. These parameters can be 
modified in LanHEP, or directly within CalcHEP once the file has been compiled. 
The values given to the parameters in Table 3.1 are the ones which were used to 
generate the cross-sections of Table 3.5 as we will explain later on. 
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Parameter 


Value 


Comment 


e 


0.31343 


Electromagnetic coupling constant (a(Mz) =1/128) 




1.21978 


Strong coupling constant {Z point) (PDG-2010) 


sin 9\Y 


0.48094 


sine of the Weinberg angle (PDG-2010) 




800 


Inverse radius of the fifth dimension 


AR 


20 


UV cutoff multiplied by R 


liR 


1 


running scale multiplied by R 


rnz 


91.1876 


mass of the Z-boson 


rrie 


0.000511 


mass of the electron 


m,. 


0.10566 


mass of the muon 


rrir 


1.77682 


mass of the tau-lepton 


rUu 


0.0025 


mass of the w-quark 


rrid 


0.00495 


mass of the d-quark 


rric 


1.27 


mass of the c-quark 




0.101 


mass of the s-quark 


^top 


172 


mass of the t-quark 


rrib 


4.67 


mass of the 6-quark 


rriH 


120 


mass of the Higgs 



Table 3.1: Parameters of the MUED at one-loop model file. Mass scales are in GeV. 



We next define the variables that are functions of the fundamental parameters of 
the model. For example, the vev of the Higgs is defined as a function of the Z-boson 
mass, the Weinberg angle and the electric charge. The first main set of 5D variables 
are the wave-function normalisation factors Zi introduced in section 2.2 and which 
are chosen to match the results of [10]. They depend on the gauge couplings -R~^, 
A and /i. From these functions, we can declare the different mixing angles found 
in section 2.2 and appendix A. For the ^("^-l^^H -W^^'^'^ -x^'''^ 

mixings, we also declare the mixing matrix entries of eq.(A.5) and (A. 9). These sets 
of variables are crucial to reconstruct the 5D fields in terms of the mass eigenstates. 

Once all the parameters are defined, we can declare the particles of the model 
(after compactification of the extra dimension) by their type (vector, spinor or scalar), 
name, mass and colour charge. Our file contains up to the second KK level. The 
reason is that, even though only the first KK level is expected to be within LHC 
reach, KK number is conserved and hence n = 2 KK particles can appear as virtual 
particles in the processes. Particles with n > 2 cannot appear in tree-level diagrams 
involving n = 0, 1 particles in the external states but would need to be considered if 
one wished to go to higher loop order or to consider higher KK modes in the external 
states. 

In Tables 3.2, 3.3 and 3.4, we list the names given to the mass eigenstates in the 

file. 
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Name 


Common name 


Particle 


Mass 


G 


G 


Gluon 




A 


7, A 


Photon 




Z 


Z 


Z boson 


MZ 


W+/W- 


W+/W- 


bosons 


MW 


el, e2, e3 


e, /i, r 


Electron, Muon, Tau 


Me, Mmu, Mtau 


nl, n2, nS 




Neutrinos 




u, d, c, s, t, b 


u, d, c, s, t, b 


Quarks 


Mu, Md, Mc, Ms, Mb, Mtop 


H 


H 


Higgs boson 


MH 



Table 3.2: The zero-mode (SM) particle names of the MUED at one-loop model file. 



Name 


Common name 


Particle 


Mass 


~G_1 




Gluon 


MGl 


~P_1 


p(i) 


P boson 


MPl 


~V_1 




V boson 


MVl 


~W+_1/~W-_1 




bosons 


MWl 


~el_l, ~e2_l 




Electron 1 & 2 


Mali, Me21 


~mul_l, ~mu2_l 




Muon 1 & 2 


Mmull, Mmu21 


~taul_l, ~tau2_l 




Tau 1 & 2 


Mtaull, Mtau21 


~nl_l, ~n2_l, ~n3_l 




Neutrinos 


MneLl, MnmLl, MntLl 


~ul_l, ~u2_l 




Up 1 & 2 


Mull, Mu21 


~dl_l, ~d2_l 




Down 1 & 2 


Mdll, Md21 


~H_1 


hw 


Higgs boson 


MHl 


~aO_l 


aO(i) 


a° scalar 


Mai 


~a+_l/~a-_l 




scalars 


Macl 



Table 3.3: First KK level particle names of the MUED at one-loop model file. The second 
and third generation of quarks are declared likewise by replacing u — )• c, t and d — )■ s, b. 
The notation with "~" for the first KK level is standard for dark matter studies. 



Having defined the 4D particles of the theory, we reconstruct the 5D interacting 
fields simply by inverting the several mixing matrices presented in section 2.2 and 
appendix A. We can then implement the Lagrangian exactly as it is written in (2.9). 
Upon compilation, the program integrates over the extra dimension and calculates 
the masses and vertices of the theory for the particles that were defined. Needless to 
say, that the Lagrangian is implemented in a 5D way ensures that gauge invariance 
is preserved through compactification. 

Let us comment on the gauge-fixing and Fadeev-Popov ghost terms in the La- 
grangian. LanHEP allows us to write a model in both unitary and Feynman-'t Hooft 
gauges by adding an option "gauge" to the declaration of the gauge bosons. The 
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Name 


Common name 


Particle 


Mass 


G_2 


G(2) 


Gluon 


MG2 


P_2 


p(2) 


P boson 


MP2 


V_2 


y(2) 


V boson 


MV2 


W+_2/W-_2 


Vr+(2)/iy-(2) 


bosons 


MW2 


el_2, e2_2 


p(2) J2) 


Electron 1 & 2 


Mel2, Me22 


mul_2, mu2_2 


(2) (2) 


Muon 1 & 2 


Mmul2, Mmu22 


taul_2, tau2_2 




Tau 1 & 2 


Mtaul2, Mtau22 


iil_2, n2_2, n3_2 


uP, uP 


Neutrinos 


MneL2, MninL2, MntL2 


ul_2, u2_2 


(2) (2) 


Up 1 & 2 


Mul2, Mu22 


dl_2, d2_2 




Down 1 & 2 


Mdl2, Md22 


H_2 




Higgs boson 


MH2 


a0_2 




a'' scalar 


Ma2 


a+_2/a-_2 


a+(2)/a-(2) 


scalars 


Mac2 



Table 3.4: Second KK level particle names of the MUED at one-loop model file. 



Goldstone boson and ghost of a given field A are implicitly defined in the progam 
as A . f and A . c respectively. In order to implement the gauge fixing and ghost La- 
grangian terms, the simplest way is to reconstruct the 5D Golstone bosons and ghosts 
of the fundamental G, B and fields out of these " . f " and " . c" . For instance, one 
may use ~P_l.f, ~V_l.f and ~aO_l in order to reconstruct B^^ (see appendix A). 
The full 5D gauge-fixed Lagrangian with Fadeev-Popov ghosts is then easily written 
in the fundamental basis. 

Finally, let us mention that we have taken particular care of the four-point 
vertices involving KK gluons. Indeed, neither LanHEP nor CalcHEP recognise the 
colour structure of these interactions. The common trick is to introduce non-physical 
auxiliary fields [34, 36]. For the interested reader, details are given in appendix B. 

3.2 Consistency checks 

At the LanHEP level, a first important validation is that of the masses: an error 
is generated if the masses declared within the model (when particles are declared) 
differ from those calculated after integration. This step ensures that the whole di- 
agonalisation procedure is correct. Furthermore, within CalcHEP, we have checked 
that our model is gauge invariant by calculating several dozen 2 — 2 and 2 — )■ 3 
processes in both unitary and Feynman-'t Hooft gauges. 

Let us now compare some cross-sections produced within the model described in 
[34] which we denote by DKM, and our implementation (BBMP). In Table 3.5, we 
list processes involving KK gauge bosons. The parameters used are given in Table 
3.1 and the center-of-mass energy is y/s = 7 TeV. Although these processes are 
unrealistic for phenomenology, they are useful to understand the differences between 
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the two models. For instance, since the Higgs sector is modified, so are the gauge 
boson eigenstates. Instead of the KK-photon 7'^"'^ and the KK Z'^"^-boson, our model 
contains the gauge bosons P'-"-' and V^"'^ introduced in section 2.2. Whenever there is 
an ambiguity concerning the initial states, we replace 7^^-' — ?■ P^^^ and Z^^^ — )■ V^^\ If 
the ambiguity is over the final states, we simply sum over all possible combinations. 



Process 


DKM 


BBMP 


g« ^GG 


3.135 X 


10 




3.135 X 10 


G ^ G 


3.183 X 


10^ 




3.183 X 10^ 




3.170 X 


10^ 




3.157 X 10^ 




3.952 X 


10-2 




3.888 X 10-2 


Gd) 7(1) ^66 


6.761 X 


10-3 




6.195 X 10-3 




6.837 X 


10-3 




5.530 X 10-3 




1.830 X 


10-2 




1.880 X 10-2 




6.910 






6.878 


w^+w ty-(i) ^ 


2.040 






2.041 


v^^+a) ty-(i) ^ Z7 


1.226 






1.226 




3.507 






3.114 


py-(i) 7 7 


1.842 X 


10-1 




1.842 X 10-1 


Z7 ^ ly-W 


1.738 






1.738 


^(1)^(1) ^ yy-(i) 


3.635 X 


102 




3.615 X 102 




1.134 X 


102 




1.097 X 102 




1.820 X 


102 




1.816 X 102 








p(i) 


5.194 X 10-2 


ly+ ly-(l) ^ ZW7 


2.858 X 10 




2.853 X 10 








tota 


2.858 X 10 








p(i) 


2.015 X 10-1 




1.134 X 10^ 




1.097 X 102 








tota 


1.099 X 102 



Table 3.5: Sample of processes with two gauge bosons for cross-section comparison (in 
pb) between previous implementation ([34], DKM) and our implementation (BBMP). 

As we can see, the two implementations are in good agreement for all processes 
considered. There are however some slight differences. Firstly, processes that involve 
both the SM and the first KK gluons give the same cross-sections, while the process 
with four C^i^ disagrees. This is due to the way the four-gluon vertices are split, as 
explained in appendix B: when a second KK level is present in the theory, the splitting 
cannot be the same as in the DKM implementation. The second set of processes 
involve KK gauge bosons into SM quark pairs. Let us focus on 7(^)7^^) tt. In 
the DKM model, electroweak symmetry breaking was not included for KK particles. 
Hence, this process is only possible via the exchange of a KK left-handed or right- 
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handed top. In our model, a (SM) Higgs can be exchanged as well, which explains 
the 25% discrepancy. The remaining processes considered in Table 3.5 involve only 
SU(2)xU(l) gauge bosons. The differences observed are again easily explained by 
the exchange of Higgs bosons H, H^^^ or a^^^\ Conversely, the processes which 
agree perfectly are those containing a SM photon simply because no Higgses can be 
exchanged at tree level. 

4 Phenomenology 

4.1 MUED parameter space, signal rates and signature under study 

Neglecting the weak dependence on the cutoff A of the theory, the parameter space of 
MUED is unambiguously defined by the inverse compactification radius and the 
Higgs boson mass mn- Moreover, this parameter space is constrained by experiment. 
Cosmological measurements of the density of dark matter limit (related to the 
mass of the lightest Kaluza Klein particle) to be below about 1.6 TeV [11]. In 
addition, tests of the Higgs sector provide a very powerful constraint on MUED. In 
cosmology, for example, the Higgs boson mass is bound from above by the requirement 
that the dark matter candidate should be neutral [11]. This limits the Higgs mass 
to be below 230 GeV. 

Constraints on the MUED parameter space also come from collider experiments. 
As we know, the discovery of a Higgs-like particle was claimed on 4th of July 2012, 
based on data released by the CMS [37] and ATLAS [38] collaborations. With 
detailed information on individual Higgs boson production and decay processes pro- 
vided by CMS and ATLAS experiments, one can understand much better the nature 
of the Higgs boson, interpret it within MUED and set limits on the model. In [19, 20], 
which established such limits, it was stressed that the signals from the Higgs boson 
from MUED are always enhanced as compared to those of SM. Therefore, the LHC 
sensitivity to SM Higgs boson would lead to even better sensitivity to the Higgs 
boson within the MUED scenario. 

One should also mention that universal extra dimensions can be limited via pre- 
cision electroweak constraints. Using results for MUED contribution to the oblique 
S,T and U parameters from [18] and recent electroweak fit of the Standard Model 
after the discovery of the Higgs-like boson at the LHC [39] {S\u=o = 0.05 ± 0.09 and 
T|[7=o = 0.08 ± 0.07 with a correlation coefficient of +0.91) in Fig. 5 we present the 
S-T parameter space allowed by experimental constraints at 95% confidence level 
(CL) (left) as well as 95% CL for the m//-i?~^ parameter space (right). 

The combined effect of all the constraints mentioned above is shown in Fig. 6 
including constraints from measurements of dark matter relic density, electroweak 
precision tests and limits from the Higgs boson search. This figure also demonstrates 
how fast constraints are evolving from the LHC side from the Higgs search: compare 
1 fb"^ (left) to 5 fb-i (right) of data. 



- 19 ~ 




Figure 5: Left: S-T parameter space allowed by experimental constraints at 95% CL 
(using S\u=o = 0.05 ± 0.09, T\u=o = 0.8 ± 0.07 and S-T correlation +0.91). Right: 95% 
CL for niH — parameter space. 
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Figure 6: Combination of Higgs constraints from 1 fb~^ (left) and 5 fb~^ (right) of LHC 
data, electroweak precision and dark matter relic density limits on the MUED parameter 
space. 

In this paper we explore the LHC potential of directly testing MUED through 
production of KK-particles. Because of the small mass split of MUED KK-particles, 
which is below about 30% of R"^ (as illustrated in Fig. 2 and Fig. 3), the main signal 
will come from strongly produced particles - KK-gluons and KK-quarks. In Fig. 7 
we present the principal decay chains of the n = 1 KK-gluon for = 800 GeV 
(top) and R-^ = 1500 GeV (bottom). 

One can see that for sufficiently large R~^ (e.g. R^^ = 1500 GeV as compared 
to R~^ = 800 GeV) the additional decay channel of G'-^^ — j- tt^i^ opens up due to 
increased mass split between C^^^ and t^i \ This channel followed by t^i^ — ta^'' 
tiu decay chain would lead to a quite striking signature of four top quarks, two leptons 
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and missing transverse momentum (remembering that G^^-* is pair produced). 

A generic feature of MUED is the high lepton multiphcity in KK quark and 
gluon decays, which can be even higher than the quark multiphcity. This can be 
seen by looking at the decay products of KK quarks and gluons in Fig. 7 which are 
dominated by orange colour, indicating the lepton decay channels. The reason for 
this pattern is in the specific mass hierarchy: 

"^JD > mw(i),mv(i) > m (i),m^(i) (4.1) 

^1,2 '■1,2 

m m > m,4/(i) yd) > m (i) (4.2) 

''1,2 ''^1,2 

which is defined by the MUED radiative corrections. 

The mass hierarchy noted in Eq. (4.1) tells us that the W^'^ and V^^^ bosons 
will always decay to leptons (SM or KK) and never to quarks. This happens because 
decays of and V^^^ to real and Z bosons are either forbidden, due to the 

insufficient mass split between y(^) and P^^\ or — when — W^P^^^ or 

are kinematically allowed (for large values R ^ ^ 1.4 TeV) — they are 
suppressed by the phase space and the small p(^)-V^(^) mixing which governs this 
decay. 

This is a very specific and unique pattern for MUED contrary to, for example, 
SUSY GUT theories. In those theories, sleptons are typically heavier than gauginos 
and therefore gauginos primarily decay to quarks because of the high SM and 
Z boson branching ratios to quarks. 

The lepton multiplicity in MUED can be easily as large as four (from the pair 
production of KK gluons or quarks). In principle, it can be as large as eight, but the 
probability of this is very low since such a process would involve rare decays of KK 
leptons into KK Higgs bosons. Moreover, leptons produced in this way will be quite 
soft because of the very small mass difference between £2^^ and a^^^^ and so they are 
unlikely to pass the experimental lepton selection cuts which we discuss below. 

It should also be stressed that the small mass split (by up to about 40% relative 
to the LKP mass) of the MUED mass spectrum, and its specific pattern, unambigu- 
ously defines the dominant production channel of MUED particles at the LHC and 
therefore dictates the search strategy: looking for evidence of MUED in the form 
of strongly produced KK quarks and gluons. The production rate of KK-gluons 
and quarks is determined by their masses which are functions of the and A 
parameters. 

Fig. 8 presents the cross sections for the G^^'^G^^\ G^^'^q^^'^ and q^^'^q^^'^ processes. 
(Here, q^^^ represents all n = 1 KK quarks, remembering that there are two KK 
quarks for each SM quark at each KK level.) One can see that q^^'^q^^^ plays a 
dominant role because it receives contributions from the multiple KK-quark fiavors. 
The G^^^g^^) production cross section is about factor 2-10 below the q^^^q^^^ one, while 
the is below the G^^^g^^) by another factor 5-20. In our analysis we sum 

over all three subprocesses. The calculation of the total cross section has been done 
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Figure 7: The principal decay channels of the n = 1 KK gluon as well as decay channels 
of its subsequent decay products for = 800 GeV (top) and = 1500 GeV (bottom). 
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Figure 8: Cross section of strongly produced KK quarks and gluons: left - for ^/s = 7 TeV 
with the decomposed contribution from G^^^G^^\ G^^^q^^^ and q^^\^^^ production; right - 
the total + G^^^g^^) + g^^^g^^) production cross section for ^/s = 7, 8 and 14 TeV 

collider energies. The CTE6L PDF was used, with QCD scale Q = A = 20R~^. 

at tree-level, so our conclusions on the LHC potential to explore MUED are on the 
conservative side since one would expect the NLO and NLL K-factors to be of the 
order of 2 and above, similar to the SUSY case studied in [40-43]. In our estimates 
we are using CTE6L parton density functions (PDF) evaluated at the QCD scale 
equal to while A was chosen to be 20R~^. 

We will show below that the LHC is sensitive to the process pp — 
(^(i)g(i) _|_ qWqW in the tri-lepton channel if it has a cross-section as large as around 
100 fb. It can be seen from Fig. 8 that this sensitivity is reached for below about 
1 TeV, 1.1 TeV and 1.8 TeV for the LHC with ^ = 7, 8 and 14 TeV respectively. 
In this paper we have focused on ^/s = 7, 8 TeV studies relevant to present LHC 
energies; however, we would like to stress that the 14 TeV LHC will significantly 
extend the limit on R~^, thus completely covering the MUED parameter space [? ] 
in combination with the Dark Matter relic density constraints. 

In Fig. 9 we present the charged lepton multiplicity distributions for signal sub- 
processes pp — )■ G^^^ G^^y G^^^ q^^y q^^'^ q^^^ and their sum for y/s = 7 (left frames) and 
-y/i = 8 (right frames) TeV before (top) and after (bottom) acceptance/selection 
cuts: 

pi > 10 GeV, \r]e\ < 2.5, Ai?,, = ^A0|. + At]^ > 0.5, (4.3) 

where i = e^,fi^ which we have summed together. For the chosen benchmark 
R^^ = 800 GeV, A = 20R~\ the mass of the KK gluon is 978 GeV while mass of the 
different KK quarks is in the range 904-925 GeV, depending on quark flavour. The 
tree-level cross sections are given in Table 4.1. 

One can see that the total cross section doubles with the collider energy increase 
from 7 to 8 TeV. 
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Figure 9: Lepton multiplicity distribution in MUED for production and subsequent decay 
of KK quarks and gluons for ^/s = 7 (left), 8 (right) TeV, before (top) and after (bottom) 
acceptance and isolation cuts. i?^^=800 GeV and A = 20i?~^. 





Cross section (fb) 


Process 




= 7 TeV 


= 8 TeV 


pp — )■ q^^'^q^^^ 


523 




1000 


pp G^i^g^i) 


205 




493 


pp G^i^G^i) 


15.8 




45.2 


Total 


744 




1538 



Table 4.1: Tree-level cross sections for pp — )■ G^'^^ G^^^ / G^^'' q^'^'^ / q^^'' q^^^ production for 
^/s = ^ and 8 TeV. 



It is worth stressing that jets and leptons are quite soft because of the small mass 
split of the KK spectrum, so one should try to keep the Pt selection cuts as moderate 
as possible. The lepton cuts can certainly be kept much softer than those for jets 
- this is related to both identification and trigger cuts. Even after quite moderate 
acceptance/selection cuts, as one can see from Fig. 9, the lepton multiplicity does 
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not go beyond 4. 

In this paper we restrict ourselves to only analysis of the tri-lepton signature 
for several reasons: 1) it is a clean signature with a comparatively low background 
level as we demonstrate below; 2) the signal rate is quite high for this signature in 
the range of that we have focused on; 3) in this paper we restrict ourselves to 
parton-level analysis which is actually is quite reliable for this signature - this is 
especially important for the background estimation. 

Although we restrict ourself to parton-level analysis, we do take into account 
realistic electromagnetic energy resolution, using a value of 0.15/^£;(GeV), which 
is typical for the ATLAS and CMS detectors, as well as their typical hadronic energy 
resolution of 0.5/ E(GeV), and perform the respective Gaussian smearing for lepton 
and quarks. 

4.2 Signal-versus-background analysis 

We consider the following backgrounds, ordered according to their relative contribu- 
tion to the tri-lepton signature under study. 

1. Di-boson production {VV, V = W^, Z) where the we took into account the off- 
shellness of one of the bosons decaying into lepton pairs; we have also included 
l^l^+jet production with > 30 GeV to take into account the leading-log 
QCD corrections and correct lepton Pt spectrum. 

2. tt production when both top-quarks decay leptonically, while the third jet could 
come from the semi-leptonic 5-meson decays. The probability of this is quite 
low (we found it to be of the order of 10~^) since the lepton from 5-meson 
decays tends not to be isolated from hadrons coming from the same decay 
because of the large momentum of the b-quark 

3. tiii and tiiuitili') production 

4. triple gauge boson production processes (l^l^V^) 

5. tiVV production 

6. background coming from 4-gauge boson production WW. 

The lepton multiplicity of backgrounds contributing to the tri-lepton signature 
is presented in Fig. 10 after applying the cuts given in (4.3). for -y/i = 7 (left) and 
8 TeV (right). One can see that indeed the VV{j) background is the dominant one. 
Its contribution to the tri-lepton signature is about one order of magnitude above the 
tt, tti£ and ttiu which in turn are bigger than the background from VW production. 
The ttVV and WW backgrounds are virtually negligible. 

In Fig. 11 we present the lepton multiplicity for the total background and the 
signal after the cuts in (4.3) for -y/i = 7 (left) and 8 TeV (right). One can see 
that after basic cuts the number of signal events is about a factor of 20 below the 
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Figure 10: Lepton multiplicity distribution for main backgrounds contributing to the 
tri-lepton signature. 

background level for the chosen benchmark of R^^ = 800 GeV and A = 20R~^: for 
^/s = 7 and 5 fb~^ one expects about 60 signal versus 1200 background events, while 
for y/s = 7 and 5 fb~^ one expects about 400 signal versus 8000 background events. 
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MUED signal vs background @LHC, Vs = 8 TeV 



background: 
signal 



4 5 6 7 8 

Number of leptons 



Figure 11: Lepton multiplicity distribution for background versus signal after the initial 
cuts (4.3), using R-^ = 800 GeV and A = 20R-^. 

Our task now is to explore kinematical variables in order to suppress background 
while leaving the signal intact. One of them is the invariant mass M^g of the most 
energetic opposite-sign leptons, presented in Fig. 12 for signal and background. In 
this and following figure we present cases for = 800 GeV and R~^ = 1200 GeV. 
For each case we present results for A = 20i?~^ and 40R~^ to give the reader an idea 
of how distributions evolve with A. One can see that increase of A slightly hardens 
the Mil spectrum since the absolute values of KK particle masses as well as the mass 
split among them increases as A increases; this was shown above in Fig. 3. One 
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Figure 12: Invariant mass of the two most energetic (highest Pt) opposite-sign leptons. 

can clearly see that the background receives a substantial contribution from the real 
Z-boson while signal does not. We thus veto M^i around niz with the cut 

\mz- Mu\> 10 GeV. (4.4) 

In Fig. 13 we present the lepton multiplicity distribution for signal and background 
after cuts (4.3) and (4.4). The background for tri-lepton signature is reduced by a 
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Figure 13: Lepton multipHcity distribution for background and signal after (4.3) and (4.4) 
cuts. 

factor of about 20 down to around the level of the signal for = 800 GeV. 

For the next step we analyse the missing transverse momentum distributions 
for signal and background which are presented in Fig. 14. One can see that the 
maximum f^j- takes place at a smaller value of for the background than for signal. 
We find that the cut 

> 50 GeV (4.5) 
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Figure 14: Missing transverse momentum f^T distribution for background versus signal 
after (4.3) and (4.4) cuts. 

is quite safe for the signal while helping to further reduce the background. 

To further increase the signal significance we use the properties of the transverse 
momentum of leptons £i, ^3 with the first-, second- and third-highest values of 
Pt respectively in each event. In Fig. 15 we plot the signal and background Pt 
distributions for each of these leptons, P^ and P^ . 

Leptons are clearly softer in signal events so the following set of the upper cuts 
on the leptons' transverse momenta help us further: 

P^' < 100 GeV; P^' < 70 GeV; P^' < 50 GeV. (4.6) 

We also make use of the effective mass M^s =^t + J2e j distribution (the sum is 
over all final-state leptons and jets) presented in Fig. 16 and apply the empirical cut 

Meff > R-^/h (4.7) 

as our final selection. The effective mass M^s is correlated with the mass of the final 
state particles, so it is harder for the MUED signal than for the background. This 
distribution is probably the most sensitive to the value of Ai? and could be further 
used for determination of the MUED parameter space in the case that we observe 
the signal. This cut is especially important for optimisation of the search for large 
values of - it allows us to suppress the background by a further factor of 2 or 3, 
leaving the signal almost intact. 

Finally, in Fig. 17 we present lepton multiplicity distributions after application 
of all the cuts - (4.3), (4.4), (4.5), (4.6) and (4.7) - together. One can see that, after 
all the cuts, the number of signal events for = 800 GeV in case of the tri-lepton 
signature is about a factor of 4 above the background for the a/s = 7 TeV case (~ 40 
signal events vs ten background events for 5 fb"^) and and about a factor of 5 above 
the background for y/s = 8 TeV (~ 150 signal events vs 30 background events for 
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Figure 15: Transverse momentum of the leptons with the first, second and the third 
highest Pt, after cuts in (4.3) and (4.4), for the background and signal. 

20 fb"^). For the R'^ = 1200 GeV corollary is at about the sinffle-event level versus 
six background events for y/s = 7 TeV, 5 fb~^ and at about the six-event level versus 
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MUED signal vs background @LHC, ^s = 7TeV 
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MUED signal vs background @LHC, Vs = 8 TeV 
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Figure 16: Effective mass distribution for signal and background after (4.3), (4.4) and 
(4.5) cuts. 

15 background events for ^/s = 8 TeV, 20 fb^^. 

We should point out that for larger (above around 1200 GeV) the sensitivity 
of number of events to variations in AR increases. This is because the cross-section 
for producing the strongly-interacting KK particles falls off increasingly rapidly with 
mass due to PDF suppression, so small changes in the mass caused by varying AR 
lead to large changes in the production cross-section. 

Another remark is in order. Looking at Fig. 17 one could get the impression 
that the di-lepton signature is a much more promising signal for MUED, but this 
is not necessarily true. Here we have explored and presented only the complete set 
of backgrounds relevant for signatures with three or more leptons. In the case of 
the di-lepton signature, one should go beyond parton-level simulations, and beyond 
backgrounds such as W/Z + jets, WW/WZ/ZZ + jets, and explore the so-called 
"QCD backgrounds". These come from, for example, W + jets production, 

where the second lepton is faked in the detector by jets (different flavours will have 
different fake probabilities) and/or photon conversion. Such study requires simu- 
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Figure 17: Lepton multiplicity distribution for background versus signal after all cuts - 
(4.3), (4.4), (4.5), (4.6) and (4.7) - for LHC @ 7 TeV (top) and LHC @ 8 TeV (bottom). 

lation probably even beyond the "fast detector simulation" level and so is outside 
the scope of the present paper. The previous study on the LHC phenomenology of 
MUED [16] has actually shown that the LHC reach for the tri-lepton signature in 
MUED is comparable to, or even better than, the di-lepton signature. We should 
stress that, in comparison with ref. [16] (the results of which we have reproduced 
quite closely for the tri-lepton signature when applying the same cuts to the same 
set of backgrounds) we suggest a more sophisticated and optimal set of selection cuts 
and consider a more complete set of backgrounds. We comment on previous studies 
in more detail at the end of the next subsection. 

4.3 LHC discovery reach 

Now we are in a position to produce results for the LHC reach. We have per- 
formed a Monte Carlo signal simulation for the grid in the {R~^,AR) plane 
using the power of the High Energy Physics Model Database (HEPMDB) 
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[http : / /hepmdb . soton . ac . uk/ ] [44] and have applied the combined set of cuts (4.3), 
(4.4), (4.5), (4.6) and (4.7) described above. 
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Figure 18: LHC @ 7 TeV (top) and LHC @ 8 TeV (bottom) exclusion and discovery 
potential for MUED for different luminosities. 




The results are shown in Fig. 18 in terms of exclusion (at 95% CL) and disovery 
(5ct) contours for ^/s = 7 and 8 TeV and different luminosities. For both criteria, 
exclusion and discovery, we define the statistical signal significance a as 

a = ^ ^ (4.8) 

VNb + Ns 

and require a > 2 for exclusion region and a > 5 for the discovery region. The 
^5(B) = <^siB)'^ denotes the number of signal (background) events for an integrated 
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luminosity L. The region on the left-hand side of the blue contour is excluded at 
95% CL, while the region on the left-hand side of the red contour is discoverable at 
the LHC. 

It is worth commenting on the shape of these contours. For the low values of 
Ai?, the LHC sensitivity contours bend to the left for low values as one can see 
by observing red discovery contours for -^i = 7 for all quoted luminosities, and also 
the red discovery contour for -y/s = 8 for the 5 fb ~^ luminosity. This behaviour is 
related to the fact that for low values of Ai? the mass split of the KK particles is 
quite small, especially for low values, as one can see in Fig. 3. The decrease of 
this mass split is correlated with the softness of the final state leptons which, in turn, 
is correlated with lower cut efficiency and respective LHC sensitivity. On the other 
hand, one can see that the LHC sensitivity contours start bending to the right in 
the larger region, even for low values of Ai?. This happens because an increase 
in leads to a large-enough mass spht to make the leptons efficiently pass the 
selection cuts. 

One can also see that, for large values of AfZ, LHC sensitivity also drops (i.e. the 
contours bend to the left). This is related to the fact that as Ai? increases the KK 
spectrum hardens which leads to a drop in the signal production rates, especially for 
large values of R~^ . 

For a total integrated luminosity of 20 fb which has been reached in the end 
of this, 2012 year at a/s = 8 TeV, the best sensitivity of the LHC is actually for 
lower values of KR. We would like stress that this is a new finding which was not 
noticed in previous studies which has generically stateed that there is a decrease of 
sensitivity in the lower KR region. 

In fact, the shape of LHC sensitivity contours - governed mainly by the values 
of the signal cross section for large values of R~^ and £ - is quite consistent with the 
shape of contours of constant mass of KK gluons and KK quarks: these contours are 
shown in Fig. 19 in the same plane of parameters. 

To conclude, we have found that the ^/s = 8 TeV LHC with 20 fb ~^ of data 
can discover MUED in the tri-lepton mode with R'^ < 1050 GeV and exclude the 
region with R'^ < 1200 GeV at 95%CL with AR in the whole range under study. 

Using this approach we estimate an optimal LHC limits for present and recent 
energies and luminosities within a consistent model of MUED. 

We would like to comment here on the results of previous studies. First of all, 
the closest phenomenological analysis to ours was performed in ref. [16]. We have 
improved it by adding additional backgrounds (most importantly Mi and Mvitiiv)) 
as well as suggesting a more optimal set of selection cuts which include the following: 
relaxation of the upper cut on the leptons' transverse momentum; application of an 
upper cut for the transverse momentum of the third lepton; and the addition of a 
final i?~^-dependent Mgff cut. These steps allow us to extend the LHC limit further 
than in ref. [16]. For example, at the 7 TeV LHC with an integrated luminosity of 
L = 2fb~^ and AR = 20, the expected lower limit on R~^ from our studies is about 
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750 GeV, versus about 700 GeV from ref. [16]. For higher energies and luminosities 
the improvement due to the cuts we have suggested is even bigger. 

We should also mention the results of ref. [17], which explored the potential of 
using the MT2 variable for the leptons+jets + I^t signature. This method definitely 
improved on the results of the previous study of the same signature [13], which used 
conventional cuts and observables. The authors of ref. [17] also demonstrated the 
advantage of their method over the four-lepton signature with conventional cuts. We 
would like to note that the tri-lepton signature that we use in our study also has 
a clear advantage over the four-lepton signature (which has a lower rate by about 
one order of magnitude). Moreover, as one can see from our results, the tri-lepton 
signature is also more powerful than the leptons -|- jets + signature with MT2 
variable in ref. [17], which predicts the LHC lower limit on to be less than 
500 GeV for L = 2 fb~^ and KR = 20 at the 7 TeV LHC. 

In the most recent study of MUED collider phenomenology [15] the authors 
successfully used so-called Rt and a shape variables, which take into account the 
topology of the of jets + signature (there are necessarily at least two quarks and 
two LKP particles in the final state at the parton level). This study has greatly 
optimised the LHC sensitivity to MUED for the jets -|- signature. For example, 
the authors showed that the 7 TeV LHC can exclude R~^ up to 700 GeV with 
L = 2fb~^. On the other hand, the overall cut efficiency is as low as 10~^, which 
leaves the tri-lepton signature as the best for setting limits on - or discovering - 
MUED. For example, using the tri-lepton signature and cuts that we suggest, for the 
same LHC energy and luminosity as above the LHC can exclude R~^ up to 850 GeV; 
this can be seen from Fig. 18. 

Finally, we mention ref. [14] which has studied the MUED signal using the top- 
quarks signature. This channel, although not being 'number one' for discovery, 
could be very useful for understanding if the properties of the underlying theory are 
consistent with MUED. 

5 Conclusions 

In this paper we have first worked out a consistent model of Minimal Universal Extra 
Dimensions and have implemented it in CalcHEP. The model is available at the 
High Energy Model DataBase (HEPMDB) in both Feynman-'t Hooft and unitary 
gauges at https://hepmdb.soton.ac.uk/ under the name "MUED-BBMP". It has 
been validated against other models for the sectors consistent with our approach [34] 
and has been demonstrated to predict different results for the Higgs sector (which 
was not complete in the previous implementations). The first part of this article, 
along with the appendices, is a self-contained review of the model and can be used 
as comprehensive documentation for the code implementation. 

In the second part of this paper we have focussed on possible signatures of 
the model at the LHC. The special MUED mass pattern (relatively compressed. 



-34- 



as compared e.g. to typical SUSY models) provides an enhancement of the lepton 
signals that can be used as a model "fingerprint" . Taking advantage of this, we have 
performed a phenomenologically-realistic study of the MUED tri-lepton signature at 
the LHC Several comments are in order. In contrast with the previous analysis, we 
have taken into account the complete set of backgrounds relevant to the signature 
under the study. We have also worked out a comprehensive set of cuts (4.3)-(4.7) 
which lead to higher signal-to-background ratio and signal significance. 

The main results of our the study are given in Figs. 18-19. In the first figure 
we show the first estimation of the LHC sensitivity to MUED using the tri-lepton 
signature for the present LHC energy and luminosity. In the second figure. Fig. 18, 
the LHC sensitivity to MUED using the tri-lepton signature is presented, which in 
particular tells us that the i/s = 8 TeV LHC with 20 fb ~^ of data can discover 
MUED in the tri-lepton mo de with R^^ < 1050 GeV and exclude the region with 
R^^ < 1200 GeV at 95% CL. 

Experimentalists have already started to analyse events with three leptons 
and missing transverse momentum, searching for production of neutralinos and 
charginos [45, 46]. We would like to encourage them to extend their analysis and 
interpret their results in the MUED framework (using our proposed strategy to opti- 
mise the cuts) to establish limits on direct production of KK quarks and KK gluons. 
It is worth mentioning that our preliminary estimations [47] show that the tri-lepton 
signature would make it possible to close the whole MUED parameter space at the 
14 TeV LHC, which is exciting news. 
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A Towards MUED at one loop: the case of U(l)xU(l) 



In this appendix, we focus on the example of two decoupled U(l)s and a Higgs 
field $ which acquires a vacuum expectation value v. We perform in detail the 
diagonalisation and gauge fixing of the Lagrangian. The starting action is gauge 
invariant but not 5D Lorentz covariant. 
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Z 



(A.l) 

where F and F' are the field-strengths of the two gauge bosons B and W respectively. 
Once again, the Z functions are intended to model the mass corrections in a gauge- 
invariant way. The covariant derivative is 
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Om-z^Bm + i^Wm^ forM = /x,5 , 



(A.2) 



where we have chosen the normalisation of the second U(l) so as to reproduce the 
SM mixing. 

The Higgs field is expanded about its vev as $ = {v + h + ix) , with h and x 
even under Z2. At the n**^ KK level the quadratic 4D Lagrangian is 
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The first line of the Lagrangian lists the kinetic terms. The second and third line 
are mass mixings. The last line is the usual derivative coupling between gauge and 
would-be Golstone bosons. 

In order for all the fields to be canonically normalised, we perform a rescaling 
5^ ^ ) and Wt^ ^ Z^'/^wt^ : 
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(A.4) 



We now proceed with the diagonalisation. 
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A.l Gauge bosons mixing 

The mass matrix of the vectors, 



Mi 

is diagonalised by 

^cos — sin 6*^^ 
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(A.5) 
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We call P^i"'' and V^""^ the mass eigenstates given by 
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with masses 
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A. 2 Scalars mixing 

The scalars mass matrix is 
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in the (^x^"\ B^\w^^'^^ basis. This matrix has a zero determinant which confirms 
that at least one of the eigenstates is a Goldstone boson. 

It is diagonalised by two rotations. We choose them to be 



23 



/l \ 

COS 923 ~ sin ^23 
\0 sin 6^23 cos 6^23 / 



with tan 2^' 



23 



2gig2/y/ZBZw 
gl/Zw - 91/Zb 



(AlO) 
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and 



O 



13 



/ COS 6^13 sin 6^13 \ 



1 



y— sin 6^13 cos6'i3y 

^(n) ^(n) 



with tan 2^ 



13 



I R2 \j Zb Zw 



We call Gr\ G^"^ and a(") these mass states. It turns out that both 2 massless 



(A.ll) 



while a^"^ has a mass 



a{n) 



(A.12) 



and the diagonal mass matrix is A^;^ = C13 C23 O23 O 



13 



/O \ 





A. 3 Derivative coupling 

In a gauge theory with Goldstone bosons, the derivative coupling has the form 
mA^d^G with m the mass of the boson and G the Goldstone boson. This coupling is 
uniquely determined by gauge invariance. The gauge fixing Lagrangian cancels out 
the derivative coupling and gives a mass mc = \/i m to the Goldstone boson in the 
-Rg-gauge. 

In our case, the derivative couplings are 
B'^-^d, [^x'-^ + ^V^^r) + (-f^x^"^ + ^^Wt^) (A.13) 

which after diagonalisation of the scalars reads 

5^ (a gS"^ + b G^"^) + W (c Gi") + d G^")) (A.14) 
with some lengthy a, b, c, d. The rotation matrix Oq 

^cos 6c — sin 6g 



/g<"'\ 



^ sin 6*0 cos 6*0 

that diagonalises the coupling (A.14) is determined by gauge invariance. We have 
introduced Gp ^ and Gq^ which are the diagonal Goldstone bosons of P;i") and qJT^ 
If we call V the initial coupling matrix, which has rank 2x3, the diagonalisation 
yields 



V 



B 



in) 



p(n) ^ .(n) 



[Oj ■ V ■ 0230130g] 



(n) 



G 
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The mass condition 




allows us to uniquely determine the rotation angle 

e„ = arctan('!^^2£§L±^) . (A.15) 

\ a cos + c sm Uy J 

where a, 6, c, d were defined in (A. 14). 

In this way, the gauge fixing Lagrangian is valid in any base, as it should. Indeed, 
the equality 

reduces toV-V'^ = It is easy to check, from (A. 13) and (A. 5) that this equality 
indeed holds. 

A. 4 Diagonalised Lagrangian 

We are now ready to proceed with the gauge fixing and write our final Lagrangian. 
The mass eigenstates are the two massive gauge bosons P*-"-* and y*-"-* ; two massless 
unphysical scalars G^p\ Gy^ which are the Goldstone bosons ; two scalars h^"'^ and 
a*-"-* . In this basis, the Lagrangian reduces to 

+mp(„)P(")9^G?^+my(„)V;(")9'^GS;) , (A.16) 



where mp(n)i "^vin) and rHain) are given in equations (A. 8) and (A. 12), and m/i(. 



n) 



ny n 



The gauge fixing Lagrangian in the gauge is meant to cancel out the last line 
of the above Lagrangian 



^GF - 2^ 



The total Lagrangian is then 



(A.17) 



-\^lin)Gf-\ml^r.)GP . (A.18) 

In the unitary gauge, ^ — >■ oo, the Goldstone bosons decouple while in the Feynman-'t 
Hooft gauge, ^ = 1, they have the same mass as the gauge bosons. 
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B Splitting of four-gluon vertices in CalcHEP/CompHEP 



In CalcHEP, it is important to take care of the vertices involving four gluons because 
of their colour structure. It was shown in [36] how to solve this issue by introducing 
Lagrange multipliers. The kinetic term F^^F^^ of SU(3) contains 



The trick is to replace it by 

^ = 



1 j 

j.a j.a iLU I „ fabcj.aiiu/^b/^c 

2^ vl ^^i^i' 



(B.l) 



(B.2) 



which reduces to (B.l) once the equations of motion for the tensor are solved and 
plugged back into (B.2). 

In the case of MUED at one-loop with two KK modes, the 5D Lagrangian is 



(5) 2 nabc nade ^{.Ld^u e , f^_G^ (5) 2 nabc nade^b Qi^id^e 

J J LL V O ^3 J J LL ^ 5 



(B.3) 



B.l Tensorial splitting 



We focus on the first term of (B.3). The expansion in KK modes and compactification 
yields 



cahc cade 
=^4 - - / 



n> 1 

^ k,l,m>l 
+ 1 ^llmnGl^'^Gl^'^G^'^^^G''^ 



k,l,m,n>l 



(B.4) 



where A^^^ and '^limn sums of Kronecker deltas arising from the integration of 
various cosines, see Appendix B in [34]. 

The function cr„KK denotes permutations over the KK index. For example 

(J ^(^^ i^) (^'^i^) (^t^ W (^'^ ^ i''^^)^ {0) Qc{k) d {I) e {m) | ^b (k) (0) ^fi d (l) e (m) 

Qb{k)Qc{l)Qfj.diO)Que{m) _^ Qb (k) qc (I) Q^J. d (m) qu e (0) 

We now truncate the tower and consider two KK levels. The Lagrangian (B.4) 
becomes 



=^4 



93 jabc jade 



n = l,2 
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It turns out that the following five auxiliary tensors are needed 



fll/ 



— no f"^" ( r7M0)(-'c{l) , nb{l)nc{0) , J_ ^6(l)^c(2) , J_ (-<fe{2)^c{l) \ 



% 

72 



71 



(B.5) 



9.r'' (^Gi«Ge(2) + -LGi(2)G:«) 



9z! 



abc 



B.2 Vectorial splitting 

We now turn to the second term in (B.3). The expansion in KK modes and com- 
pactification yield 



g2f abc fade 



J2 Gj^o^G^^^^G^-^e^^G^^") 

n> 1 

+ ^ AL^Gi(°)G^,(^^)G'^'^«G^(-) 



fc, m > 1 



k,Lm> 1 



+ 



5 ^b{k)/^c{l)^fj,d(m)^e{n 

k,l,m,n>l 



(B.6) 



where again A^^^ and A|;^^ are combinations of deltas. With two KK levels, the 
Lagrangian (B.6) becomes 



^4 = ^ glr'j""'' 



J2 (gJ/°)g^5(")g'^'^(°)g^("^ 



n = l,2 



+ G;;(l)G^(2)^M(l)G'^(2)+G'M2)^c(l)^^,(2)^e(l) 

In this case we will need vectorial auxiliary fields with the following kind of 
Lagrangian 

^ = -^^VI^V^" + Zcg^r^'V^'Gl^'^^Gl^'"^ . (B.8) 



e{2) 
5 
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We need four auxiliary vectors 



T; = ZGgzr'' [g\ 
u; = ZGgzr'' (- 




) 



(B.9) 



which, once replaced into (B.8) and eliminated, reproduce the Lagrangian (B.7). 
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